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It is shown that at the center of gravitating fluid ball there is the phase transition of second kind of 
gravitational field as a change of field's algebraic type. 
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The question about solutions of gravitational equations and researches of symmetric static 
models in General Relativity and Gravitation (GRG) today is also actual as yesterday. Petrov's 
algebraic classification of gravitational fields pQ is known as one from high-power tools of gravita- 
tional field investigations. Therefore a consideration of this problem from more general positions 
allows to analyse the fluid ball behaviour from the point of view of the Petrov algebraic classifi- 
cation. 

' In this paper is considered a spherical model of a gravitating fluid ball from more general 

base than in the papers [2J [3] . Here we suppose that our model has spherical symmetry and is 
static, i.e. all metric functions don't depend from the time variable. 

OO 

1. Metric, Tetrads and Einstein equations 

A metric interval is chosen here in Bondi's coordinates and written down as 

ds 2 = F(r)dt 2 + 2L(r)dtdr - r 2 (d6 2 + sin6 2 dip 2 ), (1) 



where functions F = F(r) and L = L(r) are the metric functions, r is a radial variable; t is a time 
variable; 8 and ip are the angle coordinates. The velocity of light and Newton's gravitational 
constant Gn are chosen equal to the unit. A determinant of a covariant metric tensor g a p, which 
corresponds to ([J) equals to det{g a p) = g = —L 2 r 4 sin 2 6. 

We shall introduce the tetrads as four orthonormal basis vectors which base on the given 
metric |T]): 

gf = y/Ffy 9(1) ^ = L5l+ l -F5l; (2) 
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5(2) M = --^( s l + isind s l); 9(3)» = --j=(Sl-isin6 (3) 
9 >{ 0) =L-H>£; gfa = 8£ + \FL-i6*; (4) 

"So = ^ + < 3) - - *>• (5) 

where i is an imaginary unit and the greek indexes run over 0,1,2,3. 

The Einstein equations in a tetrad form with a source in the form of the energy-momentum 
tensor (EMT) are written as 

G (a)(J3) = R(a)(0) ~ 29{a)(fl)R = -k7 1 ( q ) ((3 ), (6) 

where Gi a \rp\ are the tetrad components of the Einstein tensor; R( a )(p) are the tetrad compo- 
nents of the Ricci tensor; R — R , , is the scalar curvature; k = 8ir is the Einstein constant in 

' ( a ) ' 

the select system of units. The EMT of the perfect pascal fluid is written down in this case as 

T( a )(/3) = T^™^ = (fi +p)u (a) u w) - pg( a )(p) = ^u ia) u (p) +pb ia){ p), (7) 

where /i(r) is a mass-energy density; p(r) is a pressure of the perfect pascal fluid; ur a ) = 
gfocjuidx^ /ds) is the 4-velocity in the tetrad form; fr( Q )(/3) = u (a) u (p) ~ 9(a)(0) i s a 3-projector 
onto the spacelike 3-surface in the 4-dimensional world (3-metric); all functions depend from 
radial variable only. This 3-metric is orthogonal to the 4-velocity bf a \(p\u^ = 0. 

We will rewrite Einstein's equations ([6]) in a form allowing to use properties of the energy- 
momentum tensor after using a connection between the scalar curvature and a trace of the 
energy-momentum tensor R — kT : 

R (a)(P) = ~ K ( r O)(/3) - ^9(a)(0)T \ ■ (8) 

This system of equations are rewritten as the system of three equations after a substitution 
of the metric functions into the system 

£(lni)'=|( M +p); (9) 
~(lnL)' ~\(jr + l^F)' (lnF)'(lnL)') = ~xp; (10) 

->-) + |(l4)' = -4^-^ ™ 

where all derivatives are taken with respect of variable x = r/R, R is an exterior ball radius; 
X = kR 2 . 

The functions: e(x), F(x) and L(x) are connected as 



F(x) 
L{xf 



^ x ) = -rhk- (12) 



Such connection is appeared from an equality when the relational volume equals zero in a 
static case, u^.^ — 0. 

Now excluding the mass-energy density and the pressure from this system of equations, we 
get the linear differential equation with variable coefficients for function G(x): 
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G" + f(x)G' + g{x)G = 0, (13) 
where G — \/~F, f(x) — (In (p)' , <p(x) = \/e/x, and a coefficient g{x) is equal to 

I \ 2(1- 

9<P) = 2x2£ ■ (I 4 ) 
If we will input a new variable £ = C( a; ) according to a relation 

* - tSt (15) 

then the equation (|13|) goes over in the equation of a nonlinear spatial oscillator with a variable 
"frequency" fi(C(z)) 

G' c ' c + r! 2 (C(a;))G = 0. (16) 

As the equation (1161) cannot be integrated, generally speaking, in elementary functions therefore 
the "frequency" f2 can be rewritten as derivative 

«" -5(f)- 

where y = x 2 , and a function $ is an analog of Newton's gravitational potential of an interior 
part of the fluid ball. The function $ is found out of the gravitational equations through the 
function e as 

* = l-*=f / M*)^ = y v(y)Vvdy. (18) 

Besides we can easily find the relation for the pressure out of the system of gravitational 
equations (|9])- ([TTf 

XP = — % + -(l-$)(\nF)'. (19) 

x z x 

Now we will write down the function Q 2 as an expansion into series with respect to the 
variable y 

oo 

r! 2 (y) = ^a„ 2/ ", (20) 

n=Q 

here y < 1. 



Further we will choose all a„ equal to the zero, i.e. f2 2 = 0. Then (ITBl will be transformed 
into the equation 

G'U - (21) 

and function G(£) = Ci • y + C2, where C\, C2 are constants of an integration. 

In this case we have a homogeneous model of gravitational ball with a mass density /io = const 
out of the relations (fTB| and (|17p . if £7 2 = 0. The using of the connection between functions e, F 
and L (|12p and also the condition of a joint with the exterior solution of Scharzschild [3] leads 
to the well-known interior solution of Scharzschild [4]. 

Thereby we must remark that gravitational field of a spheric static distribution of the sub- 
stance belongs to the algebraic type D or the type of Petrov's classification [1] in an accordance 
with the theorem [5]. 
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In our case the gravitational field of the homogeneous distribution of a fluid has the algebraic 
type 0, i.e. the field is described by the conformally-flat solution of the Einstein equations. 

We will put n = in the expansion (f20f and now we have the constant "frequency" fl 2 = 
Qq = const. If fig is positive (Oq > 0) then a general solution of equation (fT5j) can be written 
down as a harmonic oscillating function 

G(((x)) = Goco S (n C(x) + <p ), (22) 

where ifQ is a phase displacement. 

We must here remark that the function G(C) describes spatial oscillations. 
On the other hand we find the mass density out of (fTT| and (fT5]l 

^{x)= l i Q {l-bx 2 ), (23) 
which is the parabolic distribution, and function G takes on the form 



G(x) = G cos(n (ln(-^-)/(2y/C~ 4 + y/^cj)) + <f ), (24) 

where e(x) = 1 - C 3 x 2 + C 4 a; 4 ; C 3 = xi? 2 /3; C 4 = xR 2 b/5- 

The using of gravitational equations with EMT of the perfect fluid for the given metric ([T]) 
simplifies a calculation of the Weyl tensor components. This tensor can be written in the tetrad 
components in general case as 



W / (q)(/3)( 7 )W = R (<*)(0)(>y)(6) + R h)i(a)9(l3)](6) ~ %)[(<*)S(/3)](7) _ ^ R 9(j)[( a )9(fi)](S)^ ( 25 ) 

where the square brackets denote an antisymmetric operation with the indexes. 
Mapping the Weyl tensor by 

»f =wVr feW ^" (26) 

onto the 3-dimensional euclidean space with the metric e,-j = diag(l, 1, 1), we will get the traceless 
3x3 Weyl matrix (where £kmn is an antisymmetric symbol Levi-Civita) 

W= (^)x*xdiag(2,-l,-l). (27) 

This matrix is a canonical form of Weyl's matrix of the algebraic type D according to the algebraic 
classification of spaces for the researched case and for any x ^ (n — 2m / R is a compactness, 
m is Schwarzschild's mass of the ball). 

In the point x — (at the center of ball) we have a degeneration of the Weyl matrix of D 
type into type, i.e. a continual transition x — > be accompanied by a jumplikc change of 
the algebraic type of space (D — > 0.) This change corresponds to the rank change of Weyl's 
matrix from r = 3 to r = 0. Herewith the mass density of fluid is constant and equals (io into a 
neighborhood of point x = 0. In other words the substance is not compressible nearby of the ball 
center and the algebraic type of gravitational field is the type (the conformally-flat space-time). 

Thus the obtained result is in the full correspondence with the statement of the theorem [5] . 
Moreover, the algebraic classification of gravitational fields is connected to the solution of a cubic 
characteristic equation 

A 3 + P X + q = 0, (28) 
where the parameters p and q can be found for the given traceless matrix W as 
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P=-\s P W 2 ] q = ~SpW 3 . (29) 
The equation (|28[) can be considered as the extremum condition of some " potential" function 

V(X,p, q)=V + J A 4 + X - v \ 2 + qX, (30) 

where Vq = const. 

A deformation (|30|) describes the cusp catastrophe of Whitney in accordance with the catast- 
rophe theory. Solutions of the equation (|28|) are, generally speaking, in three areas, on three 
lines and in one point of plane of the control parameters p and q when we have the static case 
(see Fig.l). 




Figure 1. The cusp catastrophe's surface and its projection onto the plane of the control param- 
eters p and q. 



It is known that the algebraic type of Riemann's space-time can be changed by the infinitesi- 
mal perturbations (herewith the rank of Weyl's matrix varies by jump) [B]. The cusp (p = q = 0) 
accords with the phase transition into type as into a most symmetric "phase" from point of 
view of the phase transitions of the second kind (Landau's phase theory [7]). Furthermore the 
parameter p is an analog of a temperature; a derivative dV /dp is the analog of the entropy; 
d 2 V /dp 2 is the analog of the thermal capacity. Last derivative varies by jump in the cusp (when 
p = ) also as the thermal capacity. 

In our case the equation ([28]) has two real coincident roots, the sum of which equals the third 
root with an opposite sign: A2 = A3 = cto = (i]/(2R 2 ))x 2 ; Ai = —2ao- The parameters p and 
q are expressed over through ao as p — 3a 2 ,; q — 2a^. A discriminant Q = (p/3) 3 + (q/2) 2 = 
describes a semi-cubic parabola (a projection of the fold lines), which corresponds to space-time 
of D type. The research case corresponds to a curve marked on Fig.l by a cross because q > 0. 

Thus a state of system has not stability in the point of an inflection of function V(X,p, q) 
with A = ao- 

The extremal values of function V(X,p, q) can be written down for all found roots of the 
equation (f28| as 

V(X 2l p) = V(X 3: p) = ^p 2 + C; V(X uP ) = ~p 2 + C, (31) 

where q and ao are expressed through the parameter p. There are the jumps of the second 
derivatives in respect of p from the function V in the point p = q = : A(V PtP ) = 1/6 and 
A(V PtP ) = —4/3. These jumps correspond to the rank jump of the Weyl matrix in the cusp [6 . 
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Hence the continual changes of the variable x (x — > 0) and of parameters p, q (p — s- 0, 
q — > 0) lead to the catastrophe: the algebraic type of space-time has the jump at the center of 
our ball model (D — > 0). 

From the viewpoint of the physical meaning there is the phase transition of second kind. The 
eigenvalues A are the parameters of the order as in Landau's phase theory [7] . 

If we will rewrite the components of the Weyl matrix through the mass density then we can 
write down Weyl's matrix ([27]) by way of 



Here one can see the transition of matrix W of the type D into the matrix of type when 
the mass density tends toward the constant value /io- 

The consideration of case with n — 1 leads to the expression ft 2 = ciq + aiy according to 
(|20|) . Then the corresponding distribution of the mass density can be written as (see (fl7 j) -(fT8 |) ) 



where c is some constant. 

The further analysis carried out similarly in a neighbourhood of the gravitating ball center 
leads to a previous result of the algebraic types' change by jump. The gravitational field has the 
algebraic type at the center (and into its small neighborhood) of the ball and in the remaining 
part of ball the gravitational field belongs to the type D. In the expansion ([2H|) there are not any 
new results for n > 1. 

2. Summary 

In the article the problem of an existence of another algebraic type at the fluid ball center than 
in the remaining part of ball is considered. The perfect fluid as the source of gravitational field 
has the arbitrary spherical distribution of the mass density. The eigenvalues task of Weyl's 
matrix is connected with the algebraic classification of Pctrov. This classification is the algebraic 
classification of gravitational field in the Riemann space-time. At the center of ball the Weyl 
matrix of spherical distribution of the mass density equals to the zero, i.e. at the center the 
algebraic type is (the conformally-flat space-time). 

On the other hand the theorem jS] says that the spherical gravitational field can be only either 
the algebraic type D or the type according to Petrov's algebraic classification. In general case 
the interior gravitational field of our fluid ball has the algebraic type D outside of the center of 
ball, where we have the type 0. 

We can see from Fig.l the cusp corresponds to the conformally-flat space-time (the algebraic 
type 0). This correspondence is strong and the algebraic type is changed by jump there. It is 
the catastrophe or the phase transition of second kind in the physical meaning. The variables 
Xi (see ([281) ) play here the role of the order parameter according to the phase theory of Landau 
[7] .And the algebraic types of the Riemann space-time play the roles of "the substance phases" 
(as in Landau's theory). 

In other words into the small neighborhood of the ball center there is the phase transition of 
the second type in the spheric gravitational field of the fluid ball for any choice of n (n > 0) in 
(|20|) as the algebraic type change of Riemann's space-time. 
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